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ABSTRACT 

We present Gemini GMOS-IFU data of eight compact low-mass early-type galaxies (ETCs) in the 
Virgo cluster. We analyse their stellar kinematics, stellar population, and present two-dimensional 
maps of these properties covering the central 5" x 7" region. We find a large variety of kinematics: 
from non- to highly-rotating objects, often associated with underlying disky isophotes revealed by deep 
images from the Next Generation Virgo Cluster Survey. In half of our objects, we find a centrally- 
concentrated younger and more metal-rich stellar population. We analyze the specific stellar angular 
momentum through the \r parameter and find six fast-rotators and two slow-rotators, one having 
a thin counter-rotating disk. We compare the local galaxy density and stellar populations of our 
objects with those of 39 more extended low-mass Virgo ETGs from the SMAKCED survey and 260 
massive (M > 10 ^°Mq) ETGs from the ATLAS^'^ sample. The compact low-mass ETGs in our 
sample are located in high density regions, often close to a massive galaxy and have, on average, older 
and more metal-rich stellar populations than less compact low-mass galaxies. We find that the stellar 
population parameters follow lines of constant velocity dispersion in the mass-size plane, smoothly 
extending the comparable trends found for massive ETGs. Our study supports a scenario where 
low-mass compact ETGs have experienced long-lived interactions with their environment, including 
ram-pressure stripping and gravitational tidal forces, that may be responsible for their compact nature. 

Subject headings: galaxies: elliptical - galaxies: dwarf - galaxies: clusters: individual (Virgo) - 
galaxies: kinematics and dynamics - galaxies: stellar content - galaxies: evolution 


1. INTRODUCTION 

Bright galaxies are known to obey well-defined sets 
of scaling relations linking their structural, kinematic 
and stellar population properties. During the past 
decade, a new generation of observational surveys 
based on optical integral-field unit (IFU) spectroscopy 
— and supplemented with spectroscopic and imag- 
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ing data covering a wide wavelength range — has 
greatly improved our understanding of these high- 
luminosity galaxies. Notable examples o f such sur- 
yeys include SAURON (de Zeeuw et al.ll2003l. ATLAS^^ 
(iCapDellari et al.l 12011 ) and GALIFA dSdnchez et al.l 

The ATLAS^d survey, which carefully examined the 
structure, kinematics, gas-content and stellar popu¬ 
lations of 260 early-type galaxies (ETGs) using the 
SAURON instrument on the 4.2m William Herschel Tele¬ 
scope, has provided many new insights into the forma¬ 
tion and evolution of galaxies belonging to this broad 
class. Employing a ATg-band magnitude selection, com¬ 
plete down of Mk < —21.5 for morphologically ETGs 
within 42 Mpc, ATLAS^^ thoroughly sampled the ETG 
population down to a stellar mass of M* ~ 6 x 10® Mq. 
As it happens, this limit coincides roughly with the 
well-known inflection point in the mass-size and mass- 
surface brightness relations that are sometimes taken 
to separate the ETGs into high- and low-m ass classes 
(see, e.g.. Figure 16 of iFerrarese et al.l 1201211 . The low- 
mass ETGs themselves appear to subdivide into sepa- 
rate branches: diffuse and compact (see, e.g., Fig. 2 of 
iDabringhausen et al.l 120081) suggesting different origins. 
Based on their simulations, diffuse low-mass ETGs are 
thought be formed mainly through tidal interactions in 
a weak field, whereas compact low-mass galaxies would 
be the low-mass counter parts of massive ETGs. It 
is difficult from imaging alone to assess the extent to 
which low-mass ETGs have properties similar to high- 
mass ETGs, and whether or not they have experienced 
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similar formation and evolutionary processes (e.g., hi¬ 
erarchical merging, violent relaxation, cold gas accre¬ 
tion, tidal stirring and/or stripping, ram pressure strip¬ 
ping, et c). Whereas strong evidence, both from obser¬ 
vations (iHuxor et al.ll2011l: ISeth et al.il2014l) and simula- 
tions (lOkazaki fc Taniguchill2000HMetz fc KrouDall2007ll 
exists about the importance of tidal stripping in the for¬ 
mation of low-mass ETGs, the respective roles of merg¬ 
ers and gas accretion is not clear yet, while the latter are 
prominent for more massive ETGs. In clusters, ram- 
pressure stripping may also be an efhcient driver for 
the transformati on of late-type star-forming galaxies into 
low-mass ETGs (IVollmer et al.ll200^ iBoselli et al.ll2009l : 

iFumagalli et al.ll2014D . 


There have been several recent attempts to provide 
an initial characterization of the kinematics and stel¬ 
lar populations of l ow-mass ETGs using long-slit or IFU 


spectroscopy (e.g., 

IMichielsen et al.l 120081: IGhilingarian 

2009; 

Koleva et al. 

2011b 

: IPaudel et al.ll2011l: IRvs et al. 

2014; 

Toloba et al. 

2014c 

). These studies have focused 


almost entirely on low-mass ETGs that have generally 
been classified as dE or dSO types (i.e., comparatively 
diffuse objects of both the nucleated and non-nucleated 
varieties). These are obvious first targets as such mor¬ 
phological classes make up the great majority of low-mass 
ETGs. 

A general finding from these studies is that low-mass 
ETGs show a surprising variety of kinematic structures: 
features that are remarkably similar to what has been 
observed in the more massive ETGs. The low-mass 
ETGs are also found to exhibit a wide range in age and 
metallicity. Thus, in terms of their stellar content, kine¬ 
matics and structural parameters, they seem to bridge 
the gap between the faintest ETGs (such as the “dwarf 
spheroidals” found in the Local Group) and the massive 
ETGs studied by ATLAS^^ . Still, the number of low- 
mass ETGs that have been examined using high-quality 
spectroscopy remains small, and existing samples include 
almost none of the compact (i.e., high-surface brightness) 
ETGs that appear to form a second branch in the pho¬ 
tometric scaling relations. Althou gh such systems ar e 
known to be very rare (see, e.g.. iBender et all I1992I) . 
a full understanding of ETGs will remain elusive until 
the nature of these objects can be understood with the 
benefit of high-quality (and preferably two-dimensional) 
spectroscopy. 

In this paper, we carry out a spectroscopic survey 
of low-mass, compact ETGs in the Virgo cluster us¬ 
ing IFU spectroscopy acquired with the 8.1m Gemini 
North Telescope. Our target galaxies have been se- 
lected from both the pioneering photographic survey of 
iBinggeli et al.l (|I985D and fro m our own Next Gene ration 
Virgo Cluster Survey lNGVS: lFerrarese et al]l20I2ll . The 
NGVS is a large (~ 900 hours) program carried out with 
the 3.6m Canada-France-Hawaii Telescope (CFHT) be¬ 
tween 2008 and 2013. The survey used the MegaCam 
instrument to perform panoramic imaging of the Virgo 
cluster, from its core to virial radius, in the u*giz fil¬ 
ters (a total area of r; 100 deg^). For a subset of the 
survey are a, r- and IL-band imaging is also available 
(see, e.g., IMunoz et al.ll20lll . The point-source com¬ 
pleteness limit for the survey is roughly g k, 25.9 mag 
(lOtr), with a corresponding surface brightness limit of 


/ig Ri 29 mag arcsec“^ (2cr above sky level). The NGVS 
image quality is also excellent, with a median j-band see¬ 
ing of FWHM Ri 0"54. A key deliverable from the NGVS 
will be an updated catalog of confirmed and probable 
Virgo cluster members, including a complete census of 
the low-mass compact ETGs that are the focus of this 
paper. 

Other papers in the NGVS series related to the gen¬ 
eral topics considered here in clude studies of the globular 
cluster populations in Virgo (|Durrell et al.ll20ll) . the dy¬ 
namical properties of star clusters, ultra -compact dwarfs 
(UGDs) and galaxies in the cluster core (|Zhu et al.ll20l4 
Zhang et al. 2015 (submitted)), a census of UGDs in 
Virgo’s three main subclusters and an examination of 
tidal stripping as a possible origin for at least some 
UGDs (Liu et al. in preparation), and investigations into 
ongoing tidal interactions of galaxies within the cluster 
envir onment (|Arrigoni Battaia et al.l[20I2t iPaudel et al.l 

[Mil . 

This paper is organized as follows. We describe our 
sample selection and data reduction procedures in 
and discuss our kinematical and stellar population anal¬ 
ysis in ^ Our results are presented in IJH and discussed 
in (jHl wherein we examine possible formation scenarios 
for low-mass ETGs in cluster environments. We summa¬ 
rize our findings in 

2. DATA 

2.1. Sample selection 

The Virgo cluster is the galaxy cluster nearest to 
the Milky Way, containing nearly 2000 cataloged mem¬ 
bers (jBinggeli et al.l 119851) . The majority of these are 
classified as low-luminosity ETGs, thus making them 
a very significant morphological group by number, and 
an important sample for understanding galaxy evolution 
in dense environmen ts. We selected our targets from 
IBinggeli et al.l (I1985D as well as fro m the NGVS. _ 

We first identified Virgo ETGs (|Binggeli et al.l[T9^ 
with Mb < —18.0, and then selected compact, high 
surface brightness objects that were objectively defined 
as lying more than Icr above the ridge-line in the 
m agnitude-surface bright ness diagram shown in Fig. 1 
of iFerrarese et al.l (j2012f) . This approach is motivated 
by the fact that more massive galaxies (M* >10® Mq) 
generally follow this relation within a Icr limit, whereas 
some low-mass ETGs lie well above the relation, a r egion 
wher e tidally-stripped galaxies may reside (e.g., iKingl 
19621: [Fdbeil 1973t iKeenan &: Innaned 119751). R e cently , 
Toloba et al. (|2014bl lcD as well as iRvs et al.l (120131) . 
have suggested that most low-mass ETGs (i.e., sys¬ 
tems with lower surface brightness than the objects tar¬ 
geted in this study) are the remains of galaxies that 
have been strongly transformed by ram-pressure strip¬ 
ping and gravitational ha rassment after falling into the 
cluster. iRvs et aP (1201411 also found that their stellar 
mass has contracted relative to their progenitors (as¬ 
sumed to be late-type). Therefore, our GFHT/Gemini 
survey — which targets some of the most compact (i.e. 
highest surface brightness) Virgo ETGs with i > 12 — 
should provide new insights into the formation of these 
extreme objects. 

We obtained Gemini data for eight of the nineteen 
galaxies that comprised our parent low-mass, compact 
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TABLE 1 

Observed and Derived Properties of CMOS Target Galaxies 


Galaxy 

(1) 

Other ID 

(2) 

RA(J2000) 

(hh:mm:ss.ss) 

(3) 

DEC(J2000) 

(dd:mm:ss.s) 

(4) 

Distance 

(Mpc) 

(5) 

d(M87) 

(Mpc) 

(6) 

d(M49) 

(Mpc) 

(7) 

mi 

(mag) 

(8) 

(mag arcsec”^) 

(9) 

Re 

(arcsec) 

(10) 

e 

(11) 

VCC 0032 

IC 0767 

12:11:02.73 

4-12:06:14.4 

16.5* 

1.394 

1.777 

13.12 

18.85 

5.60 

0.18 

VCC 1178 

NCC 4464 

12:29:21.29 

4-08:09:23.8 

15.8^ 

1.172 

0.052 

11.90 

17.91 

6.34 

0.30 

VCC 1192 

NCC 4467 

12:29:30.25 

+07:59:34.3 

I 6.32 

1.255 

0.019 

13.54 

19.10 

5.17 

0.35 

VCC 1199 

- 

12:29:34.99 

+08:03:28.8 

I 6.32 

1.236 

0.021 

14.84 

18.43 

2.09 

0.20 

VCC 1297 

NCC 4486b 

12:30:31.97 

+12:29:24.6 

I 6 . 3 I 

0.035 

1.278 

12.62 

16.72 

2.64 

0.09 

VCC 1440 

IC 0798 

12:32:33.41 

+15:24:55.4 

16.i3 

0.858 

2.092 

13.58 

20.10 

8.06 

0.01 

VCC 1475 

NCC 4515 

12:33:04.97 

+16:15:55.9 

16.7^ 

1.140 

2.420 

12.01 

18.81 

9.13 

0.35 

VCC 1627 

- 

12:35:37.25 

+12:22:55.4 

15.6^ 

0.319 

1.255 

13.89 

18.74 

3.71 

0.10 


Note. — (1) Galaxy name from the catalogue of IBinggeli et aO II1985I1 . (2) Alternative galaxy name in the NGG or IG catalogs, (3)- 
(4) J2000 coordinates (NED), (5) Redshift-independent distance when available IIMei et al.P I200'7l . ljordan et al.l ^ 120071 . IBlakeslee et al.l 
^ 1200911 or mean Virgo distance'*, (6)-(7) Projected distances in Mpc from M 87 (NGG 4486) and M49 (NGG 4472), respectively, (8)-(9) 
Apparent magnitude and surface brightness measured from NGVS i-band images, (10) Effective radius of the galaxy measured from 
the curve of growth of NGVS i-band images, (11) Ellipticity derived as a moment of the surface brightness of the NGVS i-band images, 
at i7e/2, except for two targets, VGC 1440 and VGC 1475, that we only cover up to Re /3 with our GMOS-IFU observations. 


ETG sample in the Virgo cluster, with bad weather con¬ 
ditions precluding observations of the remaining galax¬ 
ies. Despite the fact that our actual sample is not com¬ 
plete, it should nevertheless be representative and al¬ 
lows us to probe a size-mass regime that has not yet 
been surveyed with the benefit of Integral Field Spec¬ 
troscopy (IFS) data. The galaxies in our sample have sur¬ 
face brightnesses in the range 17 < < 20 mag arcsec“^ 

and effective radii of 2'.'1 ^ Re ^ 9'.'I, which corresponds 
to 165 to 740 pc a t their intrinsic di stances of, re- 
spectively, 16.3 Moc Oordan et al.ll2007l) and 16.7 Mpc 
(jBlakeslee et al.l 120091) . The properties of our sample 
galaxies are summarized in Table [TJ 

Although this was not a criterion for the selection, 
most of our objects are nucleated, or possibly nucleated, 
galaxies as shown in the ACSVCS survey (jCote et al.l 
120041200^ . They span a wide range in projected cluster¬ 
centric distance, from the center of the cluster close to 
M87 (NGG4486), out to Virgo’s outskirts (see FigH]). 
Three of them (VGC 1178, VGG 1192 & VGC 1199) have 
projected distances from M49 (NGG 4472, the bright¬ 
est galaxy of the Virgo cluster) of less than 50 kpc, 
while another (VGG 1627) is located close to the mas¬ 
sive ETG M89 (NGG 4552). A fifth galaxy, VGG 1297 
is a companion of the central massive ETG M 87 
(NGG 4486). Redshift-independent di stance est i mates 
from surface brightn ess fluctuations (|Mei et al.l I2007t 
iBlakeslee et al.l 120091) and globular cluster lum inosity 
functions (l.lordan et al.l[2?)?)7tlVillegas et al.ll2010[) . show 
that the projected distances are representative of the in¬ 
trinsic positions of the galaxies in the cluster, with the 
exception of VGG 1627, which is located at the front-side 
of the cluster at a distance of 15.6 Mpc. VGC 0032 has 
no redshift-independent distance, therefore we assume a 
mean Virgo distance of 16.5 Mpc. A more thorough dis¬ 
cussion of the environment of our sample galaxies, and 
its possible influence on their structure, kinematic prop¬ 
erties and stellar populations is presented in ^ 

2.2. Observations 

Our observations were carried out at the 8.1 m Gemini 
North Observatory, using the I FU of the Gemin i Multi 
Object Spectrograph fGMOS. iHook et '^120041) . with 
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Fig. 1.— Location of our eight program galaxies within the 
Virgo Cluster. Our targets are indicated by the large black tri¬ 
angles with green contours and are labelled by their respective 
VCC numbers. Sub-panels (a) and (b) show magnified views of 
the one square degree regions centered on M 87 (NGC 4486) and 
M49 (NG C 4472), respectively . A total of 413 “dwarf” ellipticals 
(dEs) from [iJsker et l|2007il are plotted to show the global dis¬ 
tribution of low-mass ETGs within the cluster. Colors and symbols 
correspond to their classification in sub-types (nucleated: N, non- 
nucleated: nN, disky: i and bl ue core: be). The bla ck dotted line 
indicates the virial radii, from lEerrarese et al.l II2012I1 . of the Virgo 
subclusters A (R 200 = 1-55 Mpc) and B (R 200 = 0.96 Mpc) cen¬ 
tered on M 87 and M 49, respectively. 


a combination of one- and two-slit modes due to con¬ 
figuration availability. The observations were taken in 
service mode between March 2012 and July 2013 under 
the programs GN-2012A-Q-81, GN-2013A-Q-76 and GN- 
2013A-Q-108. Table [5] provides details on the instrumen¬ 
tal setup and exposure times^ __ 

We used the GMOS-IFU (|Allington- Smith et al.ll200'^ 
primarily in its 2-slit mode, giving a 7" x 5" field-of-view 
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Fig. 2.— Composite giz-colour images from the NGVS survey (top row) and NGVS i-band residual images obtained by subtraction of 
a model computed with ELLIPSE. The ellipse fit was performed by median-sampling of the image in elliptical annuli with semi-major axis 
incremented by 10% at each iteration. During the procedure, all parameters (center, ellipticity, position angle) were allowed to vary. The 
models used to obtain the residual images shown in the lower panels do not include higher order moments and therefore only highlight 
deviations of the isophotes from pure ellipses. The green rectangle on each image indicates the GMOS-IFU FOV for the relevant mode 
(approximately b" x 1" for 2-slit mode observations, and h" x 6'' for the mosaiced 1-slit mode observations). All images are oriented to 
show north-up and east-left. The effective radius of each galaxy is indicated at the bottom right of its residual image to give an idea of the 
relative coverage of our observations. 


(FOV) with 1000 hexagonal fibers of 0".2 projected di¬ 
ameter. 500 similar additional fibers, covering an area 
of 5" X 3". 5 and pointing ^ 1 arcmin away from the sci¬ 
ence field, are used to estimate the sky background level. 
Each slit covers half of the FOV and disperses over the 
corresponding half of the detectors. To avoid spectral 
overlap on the detectors between the spectra from the 
two slits, different filters can be set in front of the grat¬ 
ing. We opted for the B600-G5307 grating combined 
with a g' filter, giving a clean wavelength range of 3980- 
5480 A. This spectral range provides extensive coverage 
of strong stellar absorption lines (e.g., H/3, Fe and Mg& 
triplet) for measuring the kinematics as well as the popu¬ 
lation parameters of age and metallicity. The instrumen¬ 
tal spectral resolution is 1688 at 4610 A for the grating 
mentioned above, with an original sampling of 0.5 A per 
pixel. We binned by a factor two in the spectral direc¬ 
tion to get a higher signal-to-noise ratio (S/N) leading to 
an actual sampling of 1A per pixel and an instrumental 
spectral dispersion of about 75 km s“^ (or a FWHM of 
2.7A at 4610A). 

Typical exposure times for our campaign were 4 x 1800 
seconds and 6 x 1800 seconds: e.g., 2 hours and 3 hours, 
respectively (see Table [5]), to reach a S/N of 20 or better 
with spatial binning while keeping enough spatial resolu¬ 
tion. One exception is VCC 1178 for which an exposure 
time of 4 X 900 seconds was considered sufficient. Due to 
scheduling constraints on the GMOS-IFU configuration, 
two targets (VGG 1297 and VGG 1627) were observed in 
GMOS-IFU I-slit mode. This mode covers only half of 
the FOV of the 2-slit mode, i.e. 5" x 3".5, and as a conse¬ 
quence the time of integration had to be doubled to cover 
the original FOV. All other parameters are similar to the 
2-slit mode except for an extended wavelength coverage, 
e.g. 4200-7200A. For consistency, however, we consid¬ 
ered only the overlapping wavelength range between the 
1-slit and 2-slit mode data in our analysis (i.e., 4200- 
5480 A). 

The GMOS focal plane consists of a mosaic of three 
2k X 4k e2vDD detectors with gaps equivalent to approx¬ 
imately 37 unbinned pixels each (^ 20 A) . To avoid losing 
spectral information, we used a spectral dithering tech¬ 
nique by shifting the central wavelength of the grating by 


TABLE 2 

Instrumental Configuration 


Galaxy 

(1) 

PAaky 

(deg) 

(2) 

PAinst 

(deg) 

(3) 

Run ID 

(4) 

Texp 

(sec) 

(5) 

VGG 0032 

76.0 

164 

1 

4x1800 

VGG 1178 

5.2 

295 

2 

4x900 

VGG 1192 

39.8 

265 

2 

6x1800 

VGG 1199 

30.8 

302 

2 

4x1800 

VGG 1297 

71.1 

1 

3 

8x2100 

VGG 1440 

81.9 

35 

2 

6x1800 

VGG 1475 

-5.8 

83 

1 

4x1800 

VGG 1627 

-81.4 

170 

3 

5x1800 


Note. _ — (1) Galaxy name from the 

IBinefgeli et al.l lUQSfiri catalogue, (2)-(3) Position 
angle of the galaxy major-axis on the sky and 
GMOS instrument in the North-East frame, (4) 

Run: 1= GN-2012A-Q-81; 2 = GN-2013A-Q-76 (2- 
slit mode); 3 = GN-2013A-Q-76 (1-slit mode), (5) 

Total on-source exposure time in seconds. 

50 A for half of the exposures, using two central wave¬ 
length settings: 4800A and 4850A. Variations in the 
point spread function between exposures, however, com¬ 
plicates the optimal combination of the spectral dithers, 
so instead we simply excluded the spectral gaps when 
performing our analysis. 

2.3. Data reduction 

We used the official Gemini IRAF package vl.11.1, re¬ 
leased in March 2012, to reduce all our data: i.e. per¬ 
form bias subtraction, flat-fielding and illumination cor¬ 
rection, wavelength calibration, cosmic ray-rejection and 
sky-subtraction. The final cubes (formed from a com¬ 
bination of individ ual exposures) wer e created using the 
XSAURON software (|Bacon et al.l 120011) . developed by the 
Centre de Recherche Astrophysique de Lyon (ORAL). We 
developed a dedicated PYRAF script to automate the end- 
to-end data reduction, ensuring a homogeneous treat¬ 
ment of the reduction procedure across all targets. 

Raw frames (e.g. science, fiat lamp and twilight expo¬ 
sures) were overscan and bias subtracted using the task 
gfprepare. The arcs were only overscan subtracted as 
they were taken in fast-read mode. The dark current of 
GMOS e2vDD GGDs is typically less than 1 electron per 
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hour and was thus ignored. 

The flat lamps within the Gemini calibration unit, 
GCAL, were used to derive the flat-field response map 
for each associated science exposure. We used the task 
gfreponse to derive the maps from the flats and we used 
one twilight per observation run (see Table [2]) to correct 
for illumination. The science spectra were extracted us¬ 
ing the aperture mask derived from their associated flat 
lamp (taken immediately before or after the exposure) 
and corrected for flat-fielding and illumination using the 
task gfextract. As explained in the § 12.21 in 2-slit mode 
each slit disperses over one half of the detector mosaic 
and a filter {g' in this case) was used to avoid overlap 
between the spectra coming from each slit. However, the 
blue end of the spectra coming from the left side of the 
FOV can still overlap with the red end of the spectra 
coming from the right side of the field. Thus, the spec¬ 
tra were shortened by 10% (^ 12 A) in the overlapping 
area to eliminate any contamination from one slit on the 
other. 

We checked the quality of the flat-fielding by apply¬ 
ing our procedure on twilight exposures, which should 
provide a homogeneously illuminated field. We found a 
residual uncertainty of about 1% between spectra from 
the same slit, but a systematic offset of 10—20 % between 
the two slits. This issue is known to arise from errors 
in the normalization from the Gemini IRAF flat-fielding 
recipe itself and not from the data. We chose not to 
correct for this offset since our analysis does not require 
accurate absolute flux calibration, and the spectral shape 
is not significantly affected. However, the isophotes de¬ 
rived from the integrated flux in the data cubes show 
some residual asymmetry as a result of this effect (see 
Fig. [5]), and so should be viewed with caution. 

We derived the wavelength solution based on the 
GCAL arc lamp exposures using the task gswavelength 
with a 4th-order Chebyshev polynomial. Due to charge 
transfer efficiency issues that can systematically affect 
the relative centroid of strong and weak lines, as well as 
introduce effects related to the read-out direction of the 
detector, we fitted only the 20 brightest lines available 
in the CuAr arc lamps. We carefully checked that those 
lines sample the whole wavelength range appropriately. 
Unfortunately, no bright lines are available redwards of 
5187 A when using a central wavelength of 4800 A which 
reduced the accuracy of our wavelength solution beyond 
that wavelength for that configuration. We therefore ex¬ 
cluded the wavelength range 5180-5480 A for the kine¬ 
matic analysis. Unfortunately, the Mg b triplet lines were 
consequently excluded from the fit for the majority of our 
galaxy sample, but the S/N of our spectra is sufficiently 
high to still allow us to reach the low velocity dispersion 
levels required (see 113.11 and Fig. [3] for details). We ob¬ 
tained typical formal RMS errors in the wavelength cal¬ 
ibration of 0.1 A and 0.08 A for the 4800 A and 4850 A 
central wavelengths settings, respectively. 

Finally, we ran the task gscrrej on each exposure to 
reject cosmic rays and cleaned the remaining bright cos¬ 
mic impacts via another cosmic ray rejection algorithm 
implemented within XSAURON. The sky contribution in 
the science spectra was subtracted by running the task 
gfskysub using the median of all dedicated sky fibers. 



Fig. 3.— Monte-Carlo simulation of GMOS-IFU observations, 
made using a model spectrum of a 12.59 Gyr, solar metallicity pop¬ 
ulation from I Vazdekis et al.l II2012I ). The spectral resolution, sam¬ 
pling and wavelength range have been matched to our observations. 
The S/N was set to 25. Gray points indicate the 2000 individual 
simulation results. Blue lines show where 68% of the points lie, as 
calculated within bins of at least 100 simulations, giving an esti¬ 
mate of the uncertainties. Errors on mean velocity, V, and velocity 
dispersion, a, are less than 11 km s“^ for velocity dispersions above 
35 km s“^. For dispersions below this value, the recovered param¬ 
eters become biased to some degree, indicated by the discrepancy 
in recovered V values. 


3. ANALYSIS 
3.1. Merging & Binning 

We used the software XSAURON to reconstruct the fi¬ 
nal cubes starting from the individual reduced exposures 
described in il2.3l All individual exposures were first 
truncated to a common wavelength range (4200-5400 A) 
and spatially aligned using isophote contours. The expo¬ 
sures were then merged by interpolating each one onto a 
common 0.2" x 0.2" grid, and then taking the optimally 
weighted sum of the exposures at each pixel. 

To increase the S/N of our data, we used 
the adaptive spati a l bin ning software developed by 
IGappellari fc Gopiiil (|2003ll . based on a Voronoi tessella¬ 
tion. The formal error spectra are not fully propagated 
through the reduction process by the Gemini IRAF pack¬ 
age. Therefore, as a first approximation, we assumed 
that the noise in a given spaxel was directly proportional 
to the square root of its integrated flux (i.e. Poissonian). 
We performed a spectral fit of the central unbinned spec¬ 
trum (see il3.2l for a description of the spectral fitting 
method used) to get a good estimate of the noise from 
the fit residuals, which we used to renormalize our Poisso¬ 
nian noise estimate before applying the Voronoi binning. 
We tested that this renormalization did not depend on 
which region of the galaxy was used to compute the fit, 
and found no dependence. 

ETGs in the luminosity range of our program galaxies 
have typical velocity dispersions of ~30 km s~^, which 
is well below our instrument dispersion of 75 km s~^. 
We used Monte-Garlo simulations with model spectra 
adapted to match the spectral properties (wavelength 
range, resolution, pixel size) of the GMOS-IFU obser¬ 
vations. Fig. [3] shows that, at velocity dispersions of 
30-40 km s“^, binning our data to a S/N threshold of 
25 was sufhcient to recover low velocity dispersions with 
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errors of ~ 11 km for both the mean velocity V and 
velocity dispersion a. This was considered an acceptable 
compromise between the kinematic errors and degree of 
spatial binning. During the binning process, we rejected 
all spectra with a S/N lower than 3 as to avoid creating 
very large poor quality bins. The adaptive spatial bin¬ 
ning on average delivers the target S/N in each spectrum 
within a standard deviation of 8%. From our simulations, 
this small spread in S/N corresponds to a variation of 
only 1-2 km s~ in the errors of velocity and velocity 
dispersion, and similarly small relative errors on the line 
strengths. 

3.2. Kinematics extraction 

To measure the stellar kinematic properties of our 
objects (mean radial velocity, V, and velocity dispersion, 
(t), we fitted the data using the penalized pixe l- fitting 
method (pPXF) from iCappellari fc EmsellemI (120041 1 
as implemented in XSAURDN. We used an empirical 
library of 89 stars tak e n fro m the MILES library 
(jSanchez-Blazauez et al.l 120061 iFalcon-Barroso et al.l 
1201111 that covers the spectral range of 3525-7500 A at 
a spectral resolution of 2.3 A FWHM (cr ~60 km s~^ 
at 5000 A), similar to (but lower than) the GMOS-IFU 
spectral resolution. We selected these stars in such 
a way to cover uniformly the parameter space {Teff, 
log(g), [Fe/H]) expected for dwarf galaxies. For each 
individual target, we first fitted the spaxel having the 
highest S/N to obtain the combination of stars that 
best reproduces the galaxy spectrum. We then used this 
single best-ht template for all spaxels to minimize the 
impact of template variations within the FOV. We fitted 
the spectra over the range 4200-5180A, excluding the 
redder wavelengths for the reasons mentioned in 1 12.31 

We performed the same kinematic extraction on the 
twilights obtained during the different runs (observed 
with the same setup and reduced in the same way as 
the science exposures) to check the quality of our wave¬ 
length calibration. Ideally, the twilight illuminates the 
full GMOS-IFU FOV and one can expect to obtain the 
same radial velocity and velocity dispersion on each of 
the 1000 science fibers. Any variations from one fiber to 
another will represent the uncertainties of our measure¬ 
ments due to, e.g., systematic instrumental signatures, 
reduction processes, etc. In 2-slit mode, we measured 
a small radial velocity systematic offset of ^5 km s“^ 
between the two slits of the GMOS-IFU instrument (i.e. 
the two sides of the field), and variations of a few km s~ 
within a single slit (for both I-slit and 2-slit modes). 
Both measurements are within the wavelength calibra¬ 
tion error (rms ~0.1 A, i.e., AU ^7 km s~^) and kine¬ 
matic errors due to noise as estimated in Fig. [S] 

We measured on the final binned maps (see § 14.11 and 
Fig. [5]) the standard deviation of the radial velocity V 
and velocity dispersion a over a 10x10 spaxels area 
where the binning was strong (i.e., at the edge of the 
GMOS-IFU FOV). We performed these measurements 
on objects having relatively flat V and/or a prohles 
(e.g., VGG 0032, VGG1440, VCC 1627) so these quan¬ 
tities represent the global true uncertainties on the kine¬ 
matics. The typical standard deviations found are of 
10 km s~^ for V and 15 km s~^ for cr, which are similar 


to the values estimated from our simulations. 

3.3. Age and metallicity measurements 

To study the stellar populations in our galaxies, we 
apply spectral fitting to derive val ues of age and metal- 
licity using the MIUSGAT models (|Vazdekis et al.l[20T^ 
that provide single age, single metallicit y stellar popu¬ 
lation models based on MILES and CaT (|Genarro et al.l 
l2001h libraries. We used a subset of 203 models cov¬ 
ering an age range of 0.5-12.5 Gyr, a metallicity range 
of —2.32 to -1-0.22 dex, and computed u sing a un i moda l 
IMF with a slope of 1.3, equivalent to a iSalnete^ (119551) 
IMF. The model spectra were retrieved from the MIUS¬ 
GAT website, sampled at 2.5 A FWHM resolution (sim¬ 
ilar to GMOS-IFU data) and covering the wavelength 
range 3464-5550A. As with the kinematics extraction, 
we used pPXF software to fit the galaxy spectra over 
the wavelength range of 4200-5400 A , excluding H/3 and 
[OIII] lines to avoid being contaminated by any possible 
emission from ionized gas. 

3.3.1. Regularized spectral fits 

As our primary approach, we used regularization of 
our spectral fits such that the best distribution of tem¬ 
plates found to fit the data is also the smoothest in the 
parameter space (in this case, age and metallicity). In 
other words, we impose the condition that the star for¬ 
mation history is the smoothest possible. The age and 
metallicity values are then derived by taking the mean 
age and metallicity of the selected templates, weighted 
by their mass contribution to the combined template. 
T his is equivalent t o the s pectral fitting ajmroach used 
by iMcDermid et 'aT] (I2015D for the ATLAS^^ survey. 

We performed this analysis on each spaxel of our 
galaxy cubes, spatially binned to a S/N of 35. This S/N 
value was chosen to preserve spatial information on each 
target while still providing robust measurements and ac¬ 
ceptable errors. Two targets, VCC 1199 and VCC 1627, 
were kept with a S/N of 25 as the data are of lower qual¬ 
ity and we opted to retain the smaller bin size. For each 
individual target, we used the spaxel with the highest 
S/N to derive the degree of regularization to use for all 
the spaxels. In this way, we keep the degree of regu¬ 
larization constant over the full FOV. The maps of the 
mass-weighted age and metallicity are shown in Fig. [5] 
and discussed in 

To estimate the errors made on the mass-weighted age 
and metalliticy, we performed the same exercise on the 
stellar population binned maps as with the kinematic 
maps and found typical standard deviations of 0.4 Gyr 
in age and 0.05 dex in metallicity. However, it is im¬ 
portant to notice that these measurements are mostly 
formal uncertainties which mostly depend on systemat- 
ics and the used spectral library. The exact impact of the 
library used for the measurements is beyond the scope of 
this paper. 

3.3.2. Best single SSP fit 

In order to compare the stellar population proper¬ 
ties of our sample galax i es with others i n the litera¬ 
ture (iKoleva et all 1201 Ibt iRvs et ahl 12014 iToloba et al.l 
l2014bUcD . we also carried out a single stellar popula¬ 
tion (SSP) analysis similar to the spectroscopic indices 
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method (|Michielsen et al.l l2007t iKoleva et al.l l2011aD . 
This method , whic h is well studied and described in 
IKoleva et al.l (j2008ll . consists of selecting the single SSP 
model that, in terms of best fits the data, taking the 
age and metallicity of this model as the measured value. 
To reduce the limitations of discreteness in the grid of 
model age and metallicity values, we interpolated the 
surface within the parameter space log(age)-[Z/H]. We 
used this interpolation to determine the position of the 
minimum and the error surface. 

We performed this analysis for each of our targets 
on the single spectrum obtained by spatially stacking 
spectra within an ellipse equivalent in area to a circu¬ 
lar aperture of Re/2, except for two targets, VCC 1440 
and VCC 1475, that we only cover up to Re/3. Effective- 
ellipses have been derived using moments of the surface 
brightness of our galaxies in the NGVS i-band images. 
The corresponding ellipticities are given in Table [T] 

We compared our res ults with published values fo r 
VCC 0032 and VCC 1192 (iSanchez-Blazauez et al.l200^ . 
VCC 1178 and VCC 147 5 (IKoleva et al.l l20lTh ) and 
VCC 1297 (|SDolaor et al.l l20I0bll bv simulating their 
long-slit data and performing the analysis on the same 
apertures. Our findings are in good agreement with these 
previous studies, as illustrated in Fig. S] The errors on 
our measurements are taken to be where the marginal¬ 
ized surface of the fit shows a difference of 1 with re¬ 
spect to the minimum value, corresponding to the 68 % 
confidence limits. 

4. RESULTS 

Several IFU-based studies of low-mass ETCs, span¬ 
ning a broad r ange of properties, have been published 
in recent years fPrugniel et al.l 120051 iChilingarianI [20091 : 
iRvs et ahlboiSH . As discussed in a our study focuses 
on the most compact low-mass ETCs, thereby probing a 
largely unexplored regime in size (165 ^ Re ^ 740 pc) 
and surface brightness /i^ < 20 mag arcsec“^. From an 
observational perspective, such objects are challenging 
targets due to their small sizes; as a result, only a few 
have had their basic kinematic and stellar population 
properties characterized. The spatial resolution and sen¬ 
sitivity of the GMOS-IFU have allowed us to obtain re¬ 
solved two-dimensional spectroscopy of the central parts 
of eight of these compact objects in the Virgo cluster. 
The kinematics, age and metallicity maps for this sam¬ 
ple are presented in Fig.[5l We also provide some further 
details pertaining to two noteworthy objects (VCC 1297 
and VCC 1475). 

4.1. Kinematics 

Six of the eight galaxies in our sample show signif¬ 
icant stellar rotation. The exceptions, VCC 0032 and 
VCC 1440, have mean stellar rotation velocities consis¬ 
tent with zero within the entire GMOS-IFU FOV: i.e., up 
to about half the effective radius for VCC 0032, and a 
third of the effective radius for VCC 1440. Interestingly, 
these two galaxies happen to reside in the outskirts of 
the cluster (see Fig. |T]). Typical mean radial velocity 
amplitudes for the other six objects are between 15 and 
30 km s~^, wh ich is similar to other low-mass, bu t more 
diffuse, ETCs (iRvs et al.ll2013l : lToloba et al.ll2011ll . Two 
galaxies, namely VCC 1178 and VCC 1297, show signif- 



Age (Gyr) - Literature [Z/H] - Literature 


Fig. 4.— Comparison of SSP a ^e and metallicity val- 
ues from our GMOS-IFU survey and IS^che z-Bla z ouez et all 
(VCC00 32, VCC 1192), IKoleva et al.l 1201 Ibl (VCC 1178, 
VCC 1475) and lSpolaor et al.ll2010bl fVCC1297'l. For the CMOS- 
IFU values, the error bars were derived from the Icr contour of 
the marginalized surface (i.e. Ax^ = 1) in the log(age)-[Z/H] 
parameter space. The upper age error of VCC 1178, VCC 1192 
and VCC 1297 are set by the boundary of the models used, 
i.e., 12.5 Gyr. 

icantly larger stellar velocities (90 and 68 kms~^, re¬ 
spectively) . Our GMOS-IFU data also confirm and high¬ 
light the presence of a counter-r otating core in VCC 1 475. 
This feature was first noted bv IKoleva et al.l (1201 IbD us¬ 
ing lon g-slit data and later re-analyzed bv iToloba et ahl 
(12014^1) (see 114.3.21 for more details). 

The galaxies in our sample have typical stellar velocity 
dispersions of 20-60 km s~^, again similar to what has 
been observed for other low-mass ETCs. Half of the sam¬ 
ple has fairly flat stellar dispersion profiles (VCC 0032, 
VCC 1192, VCC 1440 and VCC 1627). However, the re¬ 
maining four galaxies show significant central disper¬ 
sion peaks, rising to values of 60 km s“^ (VCC 1199), 
70 kms”^ (VCC 1475), 155 km s”^ (VCC 1178) and 
235 km s“^ (VCC 1297). The last three objects also show 
significant organized rotation, as well as disky isophotes 
in the photometry of the same central regions (visible in 

Fig. El) 

Fig.Elshows the composite giz-images from the NGVS 
survey as well as the i-band residual images of our target 
galaxies, obtained by subtracting from the original image 
the best ellipse fit at each radius (from the center and 
incremented by 10% at each iteration). The center, ellip- 
ticity and position angle of the fit were allowed to vary 
but the models used to obtain the residual images do not 
include higher order moments and therefore only high¬ 
light deviations of the isophotes from pure ellipses. Three 
galaxies (VCC 1199, VCC 1440 and VCC 1627) show lit¬ 
tle structure in their residual images, but typical “X- 
shape” residual of disky isophotes are apparent in the 
other five (VCC 0032, VCC 1178, VCC 1192, VCC 1297 
and VCC 1475). This is confirmed by the cosine fourth 
orde r moment in the Fourier expansion of the inten¬ 
sity (j Jedrzei ewskil 1 1 9^ which is positive out to several 
tens of arcsec for VCC 0032, VCC 1778, VCC 1297 and 
VCC 1475, and out to a few arcsec for VCC 1192. This is 
in agreement with the potential high degree of substruc¬ 
ture present in low -mass ETCs le.g.. iLisker et al.ll2007l : 

I Turner et al. II 20 i'll . Note that the lack of structure in the 
residual images can also be due to projection effects, at 
least for VCC 1440, which is a round object (e ^ 0.013) 
that does not show much rotation, and therefore could 
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Fig. 5.— Maps of mean radial velocity, V, velocity dispersion, (T, and stellar population parameters (age and metallicity [Z/H]) for 
our eight low-mass ETGs in the Virgo cluster. The box in each sub-panel indicates the values corresponding to the extreme colors of the 
colorbar on the right. The dark contours in each panel represent the iso-photometric contours of the specific target computed from the 
integrated reconstructed cube. 
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be an intrinsically flattened object seen close to face-on. 

A significant fractio n of the sample of dwarfs observed 
bv IChilingarianI (1200911 exhibited ce ntral drops in velocity 
dispersion (see also iRvs et al.ll2ni^ . In our sample, only 
VCC 1475 shows such a decreasing velocity dispersion 
profile with decreasing radius. This feature is certainly 
linked with the presence of the two counter-rotating com¬ 
ponents, which induce two peaks in the velocity disper¬ 
sion map along the major-axis, rather than a decrease in 
the center (see 94.3.21) . 

Although our small sample is not a complete census of 
the compact low-mass ETGs in the Virgo cluster, it is 
quite striking that we observe such a variety of kinematic 
signatures. Indeed, the apparent diversity is remark¬ 
ably similar to what is observed in more massive ETGs: 
from non- to highly-rotating systems, and the existence 
of large-scale, kinematically decoupled components. 

Further details regarding the dynamical support of the 
galaxies in our sample are provided in ^ via the spe- 
cif ic stellar angular mom entum parameter, A_r, defined 
in lEmsellem et al.l (j2007tl . 


4.2. Stellar populations 

Mapping the stellar population parameters of age and 
metallicity, along with kinematic substructure, can pro¬ 
vide interesting constraints on the formation and assem¬ 
bly history of galaxies. Here we present our results on 
the mass-weighted age and metallicity distributions of 
the central regions of our objects, obtained as described 
in 93. 3. II Our SSP-equivalent age and metallicity esti¬ 
mates ar e pre sented in Table [3] and discussed in more 
detail in 95.31 

On average, our objects are comprised of old to very 
old stars, with mass-weighted mean ages of 6-11.5 Gyr. 
The age variation within one galaxy spans a relatively 
large range from 0-3.5 Gyr between the center of the 
galaxy and the maximum radius covered by GMOS-IFU 
FOV, i.e. about half of the effective radius. Indeed, 
four galaxies do not show any age gradient (VCC 1192, 
VCC 1199, VCC 1440 and VCC 1627) whereas the four 
others show negative variation of 0.5 Gyr (VCC 1178), 
1 Gyr (VCC 1475), 1.5 Gyr (VCC 1297) and 3.5 Gyr 
(VCC 0032) going toward the center (i.e., younger in the 
core). Younger central stellar populat ions in low-mass 
ETG s were previously no t ed bv others (iMichielsen et al. 
2008 ; iChilingarianI 120091 : iPaudel et al.l 120111 : iRvs et al. 
201311 and are interestingly present in both non-rotating 
(e.g. VCC 0032) and rotating (e.g. VCC 1297) objects. 
In terms of metallicity, the results of the regularized full 
spectral fitting show a wide range in behavior among 
our objects, from fairly flat profiles (e.g. VCC 1178, 
VCC 1192, VCC 1199, and VCC 1440) to positive gra¬ 
dients with decreasing radius (e.g. VCC 0032, VCC 1297, 
VCC 1475 and VCC 1627). All our objects have solar 
to sub-solar stellar populations with metallicity, [Z/H], 
ranging from —0.8 to -1-0.2 dex, with a typical variation 
of, at most, 0.2 dex between the center of the galaxy 
and the edge of the GMOS-IFU FOV. In most cases, 
the central regions, which are composed of more metal- 
rich stars, correspond to yo unger stellar popu l ation s (see 
Fig. [5]), in agreement with IMichielsen et al.l (l2008f) and 
studies of more extended dwarf galaxies in the Virgo clus¬ 


ter However, VCC 1297 has an even more interesting 
stellar population structure that is discussed below. 

4.3. Specific cases 

4.3.1. VCC 1297: a central black-hole 

VCC 1297 (NGC4486b) is a compact, low-mass ETC 
that lies very close to the center of the Virgo clus¬ 
ter, just 35 kpc from M87 in projected distance. 
This M87 companion galaxy has been well studied 
and is t hought to have a central black hole of 6 x 
10® Mq (jKormendv fc Bended 119971 1. It has been re- 
porte d to show a doub le nucleus with a separation of 
0'.'15 (iLauer et al.lll996ll . but this feature was not evi¬ 
dent in imaging f rom the ACSVCS (iFerrarese et al.ll^006l : 
iCote et ^1200611 . Our observations cover the galaxy up 
to its effective radius of 2'.'3. Despite its fairly round 
isophotes (e = 0.088), a thin, fast-rotating disk-like 
structure is observed in the central 3"xl" region that 
corresponds to the location of the disky isophotes re- 
vealed by the NGVS residuals image (see Fig. and also 
IFerrarese et al.ll200CT . Its stellar radial velocity peaks at 
68 k m s~^ at ^1" along the ma jor axis, in agreement 
with iKormendv fc Bendeil (|1997[1 to within the uncer¬ 
tainties. The GMOS-IFU data also reveal that the galaxy 
rotates much more slowly (^20 km s~^) at ^ 1" away 
from the major-axis, similar to its rotation at larger radii 
(r ^ 4") along its major-axis. The velocity dispersion is 
about 150 km s~ within one effective radius but peaks, 
after a steep inward rise, at 235 km s~^ within the central 
l"xl". This strong peak in velocity dispersion leads to 
a Ser sic-dependent dynamical mass of 1.1 x 10^° Mq (see 
S I5.1l for details), implying that the black hole contributes 
more than 5% of the dynamical mass. This puts it sig¬ 
nificantly above the fundamental relation of black hole - 
host galaxy masses (|Ferrarese et al.l[2006ll . 

Very interestingly, the increase of velocity dispersion 
at the galaxy center coincides with a drop in age. The 
stars in the central regions are found to be 1.2 Gyr 
younger than the main body of the galaxy (i.e., 9.3 ver¬ 
sus 10.5 Gyr). At the same time, the metallicity map 
shows a 0.1 dex enhancement along the fast-rotating 
disk-like structure compared to the main body (0.0 ver¬ 
sus -0.1 dex). A gradient within the disk is also visible, 
with the nucleus being less metal-rich than its outer part. 
These observations suggest that an episode of star forma¬ 
tion happened at the center of this galaxy after its for¬ 
mation, possibly triggered by an external gas accretion 
that would also explain the metallicity gradient within 
the disk (i.e. more metal poor in the center). 

4.3.2. VCC 1475: a 2 — a galaxy 

VCC 1475 (NGC 4515 ) is a low-mass ETC containing a 
kinemati cally decoupled core ( KDC) that was first iden¬ 
tified by iKoleva et al.l (j2011b[ l using long-slit data and 

These values are also consistent with the well studied 
compact elliptical M32 (NGC 221) which shows a solar mean 
metallicity and a mass-weighted age of 6.8 Gyr at 2.7 ef¬ 
fective radius fl Mona^esj_e^a'n 120121 for an effective radius of 
Re ~ 100 pc llGrahamn2003r ~and an estimated dynamical mass 
of 8x 10^ Mq l|CaDDellari et aDl2006ll . Its kinematics is also sim¬ 
ilar to our GMOS-IFU objects, having a large-scale radial velocity 
amplitude of about 40 km s“^ l|Howlev et al.l[2013ll and a velocity 
dispersion ~ 60 km s“^ JCappellari et aDl2006ll . 
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Fig. 6 .— Radial stellar velocity, V, and velocity dispersion, cr, maps for VCC 1475 from GMOS-IFU d ata (top row). The x and y axes are 
both in arcseconds. The black do tted lines ind i cate the region of the galaxy observed with a long-slit bv IKoleva et al.l (I2011bl) . The bottom 
panels show comparisons of the IKoleva et al.l <l2011bl) V and a profiles with comparable long-slit data simulated from our GMOS-IFU 
observations. The error bars are computed using the light weighted standard deviation of the points considered within the simulated slit. 
IKoleva et al.l Il2011bll error bars are typically of 0.7 km s“^ and are smaller than the marker points. 


later re-analyzed bv lToloba et alJ (l2014all . This galaxy is 
located in the outskirts of the cluster, about 1.1 Mpc in 
projected distance from M87 (i.e., 0.7 the subcluster A 
virial radius). It has a mildly flattened shape (e ~ 0.2) 
within its effective radius of 9'.'3 but becomes increas¬ 
ingly flattened (e ~ 0.42) within its central 2"x2", re¬ 
sulting in an ellipticity of 0.35 within the GMOS-IFU 
FOV. The isophotes in that region are disky, as empha¬ 
sised in the residual image presented in Fig. [2] Our two 
dimensional data cover the full KDC and reveal a true 
counter-rotating disk (aligned with the major-axis) that 
has a minimal size of ^ 5" (i.e. ^ 400 pc) in diame¬ 
ter. The two stellar counter-rotating components have a 
mean velocity amplitudes of 8 km s~^ (inner component) 
and 16 km s~ (outer component, within the GMOS- 
IFU FOV). The stellar velocity dispersion map shows 
two off-centered peaks along t he major-axis, makin g 
VCC 1475 a bonafide 2—a galaxy (iKrainovic et al.l[2011h . 
The amplitudes of the peaks are of order 100 km s“^ and 
correspond to the apparent transition regions between 
the two counter-rotating disks. The velocity dispersion 
drops to 70 km s~^ in the nucleus. We thus confirm 
the dis covery of a KDC in VCC 1475 by IKoleva et al.l 
(j2011bfl , but clearly reveal its counter-rotating and 2 — cr 
nature with our GMOS-IFU data. Figure [6] presents a 
magnified view of the kinematic maps for VCC 1475 as 
well as a comparison of the rotation curve of iKoleva et al.l 
(|2011bll with that extracted from our GMOS-IFU data. 


Like VCC 1297, this 2 — cr galaxy has a stellar popula¬ 
tion that is 1 Gyr younger in its core region than in the 
external region covered by our observations (9.5 versus 
10.5 Gyr). This outward increase in age is accompanied 
by a decline in metallicity going from the central to outer 
regions (-0.1 versus -0.5 dex, respectively). These differ¬ 
ences in stellar populations seem to spatially match the 
two kinematically distinct components and are very sim¬ 
ilar to the behavior seen in VCC 0032, both in terms of 
age and metallicity gradients. The kinematically distinct 
inner component in VCC 1475, coupled to its younger 
stellar populations, strongly suggests that a dissipative 
accretion event has occurred in this galaxy, followed by 
an episode of star formation. 

5. DISCUSSION 

In §[3]and §|4l we have analyzed two-dimensional spec¬ 
troscopic data for eight compact dwarf ellipticals in the 
Virgo cluster obtained with the GMOS-IFU instrument. 
In this section, we discuss our findings in the context of 
possible ETC evolution scenarios. 

5.1. Mass, size and local galaxy density 

A key goal of our study is to understand whether the 
compactness of our program objects is related in any 
way to a specific physical evolutionary process. Many 
previous investi gations have emphasized the importanc e 
of harassment (jMoore et al.l 119981 : iBoselli et al.l 1200811 , 
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strangulat ion (ISmith et al.|| 2012 f |, ram pressure and tidal 
stripping (iGnedinI 200,11 Maver et al.l 120061 : iLisker et al.l 
120131 : iToloba et al.ll2014ch and. more generally, interac¬ 


tions with the underlying galaxy population or inter- 
galactic medium, i n the evolution of low-mass galaxies 
(jBoselli et al.l[201^ . In that spirit, we first turn our at¬ 
tention to the relationship between the basic properties 
of our targets and the galactic density in the regions in 
which they reside. 

In Fig. 13 we use the unique coverage and depth of 
the NGVS dataset to examine this relationship for Virgo 
galaxies using a mass-size diagram (i.e., effective radius, 
i?e, versus dynamical mass, M). We compute masses for 
our objects using the virial theorem, 


M = 0.1 Pin) <7R/ G 


( 1 ) 


where I3{n) depends on the Sersic index, i.e., the radial 
pr( 3 file of the surfac e brigh tness , following equatio n ( 20 ) 
of iCappellari et all (1200^ (cf. iBertin et al.l l200^ . aR 
denotes the stellar velocity dispersion within an aper¬ 
ture of radius R and G is the gravitational constant. 
We used the Sersic indices measured from the i-band 
NGVS images and derived the velocity dispersion within 
the largest ellipse covered by our ob servations to es¬ 
timate a-Rg using the equation ( 1 ) of ICappellari et al.l 
(1200^ . For an aperture of one effective radius, this 
mass is representative of the enclosed dynamical mass 
(ICappellari et al.ll2006L l^llall . Our derived masses are 
presented in Tabled 

We then combined our sa mple with 39 ( on ave rage) 
more extended objects from iToloba et al.l (l2014bll ob¬ 
served with long-slit part of the SMACKED survey, and 
representative of the low-mass ETGs of the Virgo clus¬ 
ter within the magnitude range —19.0 < Mr < —16.0. 
Dynamical masses f rom these objects were derived by 
IToloba et al.l ()2014bH in a similar way to ours, using the 
Sersic-dependent virial estimate. We also compared the 
NGVS effective radius values with those of iToloba et al.l 
()2014bfl for the objects in common and found a very good 
agreement. Therefore we can be co nfident in the relative 
position of our sample compare to IToloba et al.l (|2ni4br i 
sample in the mass-size diagram. We also extended the 
analysis to significantly higher masses using 56 Virgo 
ETGs observed as part of the ATLAS^^ project. The 
sizes and masses for the ATLAS^^ sample come from 
ICappellari et al.l (l2013al lbh . The masses were derived 
via anisotropic Jeans dynamical modeling, which were 
shown to give good a.greement with Sersic-d ependent 
virial mass estimates (ICappellari et al.l l2013afl . As ex¬ 
pected, Fig. [7] (panel a) shows that our sample contains 
low-mas s ETGs that are mo re compact than those tar¬ 
geted in IToloba et al.l[2ni4bl . Indeed, for the same dy¬ 
namical mass ( 10 ®- 10 ^° Mq ) our galaxies have effective 
radii that are nearly an order of magnitude smaller than 
their extended counterparts: thus, our study probes a 
new parameter space for low-mass ETGs. Clearly, it will 
be of interest in future IFU surveys to target galaxies 
in this under-represented part of paramete r space. Note 
that th e prominent rap in mass between IToloba et al.l 
(12m 4H1 and ATLAS® ETGs in Fig. [7] is purely due to 
sample selection effects. Two objects from our sample, 
namely VCC 1178 and VCC 1297, are of particular inter¬ 
est as they overlap in mass with the lower-mass range 
of the ATLAS^^ objects (see also Fig. (TOl), but are lo¬ 


cated within the ‘zone of exclusion’ (ZOE) defined by 
(jCappellari et al.l[2ni3b[l . Both objects have an absolute 
K-band magnitude {Mk=-21.A and -20.8, respectively) 
that is fainter than the ATLAS^'^ sample selection def¬ 
inition {Mr < —21.5). At face value, this implies an 
unusually high mass-to-light ratio for these objects, but 
detailed modeling would be required to determine the 
cause of this, which is outside the scope of this paper. 

In Fig. [3 the galaxy density is expressed in two dif¬ 
ferent flavors: ( 1 ) the luminosity of the brightest galaxy 
among its 15 closest neighbors, normalized by the lu¬ 
minosity of the target galaxy, Lmax/L] and ( 2 ) the 
surrounding galaxy density, E 15 , defined as the num¬ 
ber of galaxy per square degree within a projected re¬ 
gion that includes the 15 closest neighbors to the target 
galaxy. These two values for our sample galaxies are 
presented in Table [S] The first result that emerges from 
these plots (Fig. [3 is that the five low-mass galaxies 
that depart most dramatically from the main mass-size 
sequence (VCC 1178, VCC 1192, VCC 1199, VCC 1297, 
VCC 1627) are all found in regions of high galaxy den¬ 
sity (i.e., high E 15 , see panel b) and, most importantly, 
that also contains a massive galaxy (i.e., high Lmax/L, 
see panel a), as expected since they are in close prox¬ 
imity to either M87 (NGC4486), M49 (NGC4472) or 
M89 (NGC4552). The fact that most galaxies above 
10 ^° Mq have low Lmax/L parameters is a natural con¬ 
sequence of the galaxy luminosity function: i.e., the prob¬ 
ability that such a galaxy will be found next to an even 
more massive galaxy is lower than for low-luminosity ob¬ 
jects. In the low mass regime (< 10^° Mq), we indeed see 
a much wider range of Lmax!L values, which makes the 
high Lmax!L value of the more compact low-luminosity 
galaxies even more relevant. Interestingly, these five ob¬ 
jects have very small nu mbers or no globular cluster sys¬ 
tems (jPeng et al.ll2008D . suggesting indirectly that they 
could have experienced tidal stripping events. If, rather 
than using LmaxIL, we use the normalized total lumi¬ 
nosity Ltot/L, the results do not change. 

To better reveal the first order relations in the three- 
dimensional space, we also smoothed the data using the 
two-dimensional Locally Weighted Regression (LOESS) 
method ( Cleveland fc Devlin 1988) , implemented as de¬ 
scribed in ICappellari et aP (l2013bl) . Results are shown 
in the two bottom panels of Fig. [3 (panels c and d). This 
analysis used a regularization factor of 0.3. 

Ultimately, we performed the same exercise with a var¬ 
ious number of closest neighbors N (i.e. N=3, 10 and 25) 
and we did not see any significant changes in the trends 
except a systematic increase of LmaxIL with N (as ex¬ 
pected from the cluster luminosity function). 

5.2. Rotational support: the Xr — e diagnostic 

lEmsellem et al.l (|2011l l used the specific stellar angular 
momentum versus ellipticity diagram (Xr — e) to clas¬ 
sify massive ETGs, reporting an empirical criterion with 
which to identify slow and fast rotators (that represent 
^ 15% and 85% of the ATLAS^° sample, respectively). 
This scheme has been used subsequently to constrain 
the properties of possible low-mass ETGs progenitors 
(|Rvs et al.ll2013l : IToloba et al.l l2014cD. Accordingly, 
followed the methodology given in lEmsellem et al.l (120071) 
to derive Xr profiles for all galaxies in our sample (up to 
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Fig. 7.— Location of Virgo galaxies from our GMOS-IFU sample, those of IToloba et al.l Il2014bll observed with long-slit, and 
ICappellari et aO Il2013bl l from the ATLAS®^ survey, in the mass-size (M — Re) plane, with color-coded symbols representing a proxy 
of the local galaxy density: the a and c (left) panels use the luminosity of the brightest galaxy among the 15 closest neighbors of the target 
galaxy, normalized by its luminosity, L-maxIL ; the ^ and d (right) panels use the S 15 parameter defined as the number of galaxies per 
square degree within a region that includes the 15 closest neighbors to the target galaxy. The bottom panels show a LOESS version of the 
original dataset presented in the top panels using a regularization factor of 0.3. The gap between the low- and high-mass ETGs is due to 
sample selection effects. 


the maximum radius coverage of the GMOS-IFU FOV). 
Our FOV limit corresponds to i?e/2 for all galaxies in 
our sample except VGG 1440 and VCC 1475, which are 
only covered up to i?e/3. The measured values of e are 
presented in Table [H and Xr in Table. |3l The Xr — e 
diagram is shown in Fig. |8l 

We find six fast-rotators (75%) among our sample 
of eight galaxies: VCC 1178, VCC 1192, VCC 1199, 


VCC 1297, VCC 1440 and VCC 1627. The two remain¬ 
ing galaxies are classified as slow-rotators: VCC 0032 
and VCC 1475. The fraction of slow-rotators in low- 
mass compact ETGs (25%) is slightly higher than that 
found among high-mass ETGs, but given the small num¬ 
ber of objects in our sample, the difference is not sig¬ 
nificant. Still, it is worth pointing out that the fraction 
of slow-rotators among the entire low-mass ETGs popu- 
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TABLE 3 

Derived Kinematic Parameters, Ages and Abundances 


Galaxy 

i?e 

n 

Xr 

O-R 


Me 

Age (SSP) 

[Z/H] (SSP) 

Sl5 

Rmax / R 


(pc) 



(km s“l) 

(km s 1 ) 

(M©) 

(Gyr) 

(dex) 

(gal.deg®) 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

VCC 0032 

447.9 

2.39 

O.O 915 I 

30.9 

29.5±0.72 

(4.5±0.22) X 10* 


-0.687l°;®f3 

13.43 

3.04 

VCC 1178 

486.0 

2.52 

O. 3959 I 

138.3 

132.1±3.21 

(9.6±0.46) X 10® 


-0.052lO-°®i 

82.99 

61.89 

VCC 1192 

408.9 

1.81 

0.33581 

56.3 

53.7±1.31 

{1.4±0.07) X 10® 

10.41®;® 

0.0381®;)^® 

98.54 

283.1 

VCC 1199 

165.3 

1.19 

O. 2657 I 

53.4 

51.1±1.24 

(5.5±0.27) X 10* 

Q 1 +3-4 
^•^-3.6 

0-22l°;V 

93.41 

968.2 

VCC 1297 

208.6 

2.45 

O.I 599 I 

228.1 

217.9±5.30 

(1.1±0.05) X lOi® 

-*--*-•^-1.8 

U.zz_Q_i82 

133.50 

122.35 

VCC 1440 

629.1 

3.14 

0.1672^ 

34.8 

32.4±1.25 

(7.0±0.53) X 10* 


-0 415+0-096 

U.4iO_Q 

30.29 

3.24 

VCC 1475 

739.5 

3.32 

0.0657^ 

82.8 

77.0±2.96 

(4.5±0.35) X 10® 

12.5in 

-0.415ini7 

25.64 

0.11 

VCC 1627 

280.7 

1.62 

O. 3155 I 

40.5 

38.7±0.94 

(5.2±0.25) X 10* 

6.5l|J. 

0.038l®;i«i 

99.28 

85.82 


Note. — (1) Galaxy name from [^nggeli et al.l II1985I) catalogue, (2) Effective radius in pc derived from the effective radius and 
distance indicated in Tabled (columns 10 and 4, respectively), (3) Srsic index from NGVS Aband images, (4)-(5) Specific stellar angular 
momentum \r and velocity dispersio n within an aperture o f Rel‘2^ or (6)-(7) Estimated velocity dispersion and dynamical 

masses at at one effective radius using fUappellari et al.l II2006I) . (8)-(9) SSP-equivalent age and metallicity measured within an aperture 
of Rel‘2^ or RelS^. The errors on our measurements are taken to be where the marginalized surface of the fit shows a difference 
of 1 with respect to the minimum value in the log(age)-[Z/H] parameter space. (10) Galaxy surface density S 15 , i.e., within a radius 
including the 15 brightest neighbors, derived from NGVS data, (11) Luminosity ratio of the brightest galaxy among the 15 neighbors 
of the galaxy and the galaxy. 


lation (i.e., the iRvs et al.ll201^ iToloba et al.ll2014d and 
present GMOS-IFU sample combined) is still at least 
25%, or 10% higher than in high-mass ETGs. Inter¬ 
estingly, VCC 0032 shows gradual but large variations in 
both ellipticity and position angle within its inner 40", as 
well as disky-isophotes. It may therefore represent a true 
slow rotator, as opposed to a more face-on fast rotator. 
The low Xu value of VCC 1475 is not unexpected since 
we have shown it to be an unambiguous 2 — a galaxy 
(|Krainovic et al.ir201lL see S 14.3.21) . 

Our sample galaxies are found to have 0.1 ;< ^Re/'i ~ 
0.4 and relatively low ellipticities: e < 0.4. This is 
similar to what has been found for m ore e xtended low- 
mass E TGs by both iRvs et aP (1201411 an d iToloba et al.l 
l|2014cn . Note that the IRvs et al.l (I20T^ Xr values are 
derived at i?e, creating a systematic, but not significant, 
positive offset. None of our targets reach the high values 
of ^R ^/2 > 0.4 or e > 0.4 that are commonly observed in 
h igher mass ETG sam ples (jEmsellem et al.ll2011[) . 

IToloba et al.l l|2014cr ) suggested the existence of a cor¬ 
relation between Xr and the projected distance to the 
center of the cluster (e.g., M87) as well as the pres¬ 
ence of subtle disky structures in fast-rotators. Among 
our targets, VCC 0032, VCC 1475 and VCC 1440, are the 
outermost galaxies. They are the most extended objects 
and are located in regions of relatively low galaxy den¬ 
sity. The five remaining objects (VCC 1178, VCC 1192, 
VCC 1199, VCC 1297 and VCC 1627) are found close to, 
either M87 or a large companion (see >i5.1l) and are all 
fast-rotators. More interestingly, they are faster rota¬ 
tors than the three remaining objects and are also the 
most compact objects of our sample. While these results 
may appear in contr adiction with the trend suggested by 
IToloba et a,l.l (120143) — tha.t, on average, Xr. normaliz ed 
bv a/c Isee lEmsellem et al.ll2011l: IToloba et ani2014(Il is 
lower at small distances from M87 — their trend seem 
weaker, or even absent (see their Figure 5), for galax¬ 
ies with disk-like structures. Three of the five compact 
galaxies may have such disk features, suggested by their 
disky isophotes, so it may not be surprising if they do 



Fig. 8. — \r — e diagram. Values have been derived at i?e/2 
for all our GMOS-IFU objects except VCC 1440 and VCC 1475, 
derived at Re/3. Xjie. values of m o re ext ended low-mass ETGs in 
the Virgo cluster from [Ry;^^aL|_J20^n are sh own as green cir¬ 
cles and Xjie /2 values from IToloba et aLl Il2014d^ as blue triangles. 
For lToloba et al.l l|2Q14dj and IRvs et al.l 1)20141) samples, open sym- 
bols represent obje cts with underlying disky-structure identified in 
ILisker et al.l 1120071) . For GMOS-IFU sample, open symbols rep¬ 
resent objects with disky isophotes identified in the NGVS. The 
magenta line represents the empirical relation between Xr and e 
for an ellipsoidal object seen edge-on and having an anisotropy pa¬ 
rameter of = 0.65 X e. The black dashed lines correspond to this 
same relationship for an object seen edge-on to face-on. The orange 
line (respectively the dashed orange line) represents the empirical 
relation Xr = 0.265 x yj^ej^ (respectively Xr = 0.31 x y^) de¬ 
fined bv lEmsellem et'a n i20TiD which separates slow-rotators from 
fast-rotators. 
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Fig. 9. — Location of our GMOS-IFU galaxies (circle symbols) in the mass-size diagram combined with those of IToloba et al.l Il2014bll 
(triangles) observed with long-slit. The color of each symbol corresponds to the SSP-equivalent age (panel a) and metallicity (panel b) of 
the combined spectra within Rej'l, as indicated by the color bar in each panel. 


not participate in such a trend. Observations of larger 
galaxy samples are definitively needed to confirm, or not, 
these possible trends. 

5.3. Stellar content and star formation history 

We now analyze the age and metallicity of the stel¬ 
lar populations in our target galaxies, as tracers of 
their formation and assembly history. Because the SSP- 
equivalent ages and metallicities were derived from the 
stacked spectra, they reflect the single stellar popula¬ 
tion that dominates the light within the central regions, 
i.e i?e/2 or Re/3 (see Table [S]). No trend is found be¬ 
tween the stellar population parameters and the loca¬ 
tion of our targets within the cluster (the local galaxy 
density) or thei r Ar parameter. Thi s is qualitatively in 
agreement with iToloba et ahl (l2014bD . It is, however, in¬ 
teresting to note that all objects with underlying disky 
isophotes, aside from VCC0032, show the oldest single 
stellar population. 

As emphasized in 1)5.11 the galaxies in our sample 
are more compact than those examined in other stud¬ 
ies of low-mass ETGs. Thus, we investigated the po¬ 
tential link between their compactness and their stel¬ 
lar population age and metallicity as shown in Fig. O 
Here we present the GMOS-IFU sample in terms of their 
SSP-equivalent age and metallicity distributions in the 
mass-size diagram. To compare our sample to more ex- 
tended low-mass FT Gs, we have added the objects from 
iToloba et al.l (|2014bD . shown as triangles. No clear trend 
appears between the age and the location of the objects 
in the mass-size plane although, on average, the com¬ 
pact FTGs are older than less compact systems (panel a 
of Fig. O. A possible trend is visible between the metal¬ 
licity and the compactness of low-mass FTGs (panel 
b of Fig. [9]): e.g., more compact objects tend to be 
more metal-rich. In the context of low-mass ETGs be¬ 


ing for med through stripping and harassment. iB.vs et al.l 
(j20I4l) . showed that these objects become more com¬ 
pact than their progenitors, which are often assumed 
to be late-type, low-mass spiral galaxies. The GMOS- 
IFU sample includes mor e compact objects than those in 
iToloba et al.l (l20I4cl) and iRvs et all (12014D , so we would 
expect our objects to have been further processed. Our 
findings of older and more metal-rich stellar populations 
in the most compact ETGs is fully consistent with this 
picture. 

In the broader context of galaxy formation, there 
exists a longstanding debate on whether low- and 
high-mass ETGs show a dichotomy in their proper¬ 
ties, i.e., whether they represe nt two distinct pop- 
ulations (as first advocated by iKormend^ I1985D . or 
whether they form a more continuous sequence in 
terms of their structure, kinematics and stellar popu¬ 
lations. Most recent studies of photometric scaling rela¬ 
tions and isophotal parameters sugge st a continuity be- 


19971 IGraham & GuzmanI 12003 

: IGavazzi et al.l 

20051 

Ferrarese et al.l 120061 Gote et al.l 

2007fl. but an alterna- 

tive view is favoured bv [Kormenc 

v et al.l (|2009[) . Recent 


surveys, especially those using IFU spectroscopy, have 
characterized the detailed kinematics and stellar popu¬ 
lation properties of more massive ETGs. Using them 
as a reference sample, we have combined the 260 mas- 
sive ETGs (IO^^Mm < M < 10^^ M^) from AT LAS^° 
(|Cappellari et al.l l2013al iMcDermid et al. 120151). with 


the 39 low-mass ETGs from IToloba et ah (l2014bD and 


our eight compact low-mass ETGs. We focus here on the 
SSP-equivalent stellar population information (e.g. age 
and metallicity) as projected onto the mass-size diagram 
(see Fig. [IQI). As in Fig. [3 we smoothed the raw data 
(panels a and b of Fig. [TUI) using the LOESS method 
with a regularization factor of 0.3 to reveal the first or- 
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Fig. 10. — Location of our GMOS-IFU galaxies (circle symbols) in the mass-size diagram. Samples from IToloba et al] Il2014bl l observed 
with long-slit (triangle symbols), and lMcDermid et al.l Il2nifil ~l from the ATLAS*® survey (square symbols), had been added to the diagram. 
The colors of each symbol correspond to the SSP-equivalent age (left panel) and metallicity (right panel) as indicated by the color bars 
in each panel, respectively. The light-blue dotted lines represent the lines of constant velocity dispersion cr=20, 30, 50, 100, 200 250 

km s“^ (from left to right) as implied by the virial mass estimator M = bRecr^IG. The gap between the low- and high-mass ETGs is due 
to sample selection effects. 


der relationship in the three-dimensional space (panels c 
and d of Fig. fTOI) . 

The first striking result from this figure is that there 
is a level of continuity in age and metallicity for ETGs 
between 10 * M q and 10^^ Mq . As observed in 
iMcDermid et all ()201,5[ ). at fixed size, more massive 
galaxies are older; alternatively, at fixed mass, larger 
galaxies are younger. In term of metallicity, low-mass 
ETGs that are more compact are also more metal- 


rich, and more massive galaxies in general tend to have 
more metal-rich stella r population ([Thomas et alJl^OlCt 
IMcDermid et alJl2015f) . It is also interesting to note that 
the ages and metallicities of the most compact, low-mass 
ETGs are more representative of more massive ETGs, 
rather than of less compact galaxies with comparable 
mass. This similarity between compact low-mass and 
massive ETGs has already been observed, in term of 
metallicity only, by iChilingarian et ^ 120091 for a sam- 
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pie of 21 compa ct ellipticals in the loca l group. This is in 
agreement with lMcDermid et al.l (1201511 who showed that 
the stellar population parameters of massive ETGs in the 
mass-size pla ne roughly follow lines of iso-velocity disper¬ 
sion (see also lCaPDellari et “3111201 111 , as indicated by the 
light-blue dotted lines in Fig. [TO] In fact, we find that 
this result can be extended to low-mass ETGs as shown 
in Fig. [10] (albeit with a large scatter in the age distri¬ 
bution), with the transition between low- and high-mass 
ETGs being smooth. Taken at face value, this naively 
suggests that low-mass ETGs (10® < M < 10^° M©) 
may have been shaped by the same secular and/or en¬ 
vironmental mechanisms as their massive counterparts 
(M > 10^° Mq). However, it does not necessarily imply 
that galaxies evolve along iso-velocity dispersion lines, 
neither constant metallicity lines. Only detailed simula¬ 
tions focused on stellar population and kinematic param¬ 
eters would give us clues about the evolutionary track in 
each scenario (e.g. tidal stripping, ram pressure strip¬ 
ping, harassment, etc.). Fig. 1101 mostly confirms that 
velocity dispersion is a good predictor of the age and 
meta l licity of galaxies (Graves fc Fabeill2010l : I Wake et'H] 


(Ml IMcDermid et al.^oTlT even for low-mass ETGs. 

Finally, if we consider the two-dimensional, mass- 
weighted age and met allicity maps for our objects 
(see Fig. [5] and § 14.2L we see that nearly half our 
program objects contain a younger and more metal- 
rich stellar population in their cores. This result 
is in agreement with what has b een previous l y ob¬ 
served in “dwyF’ galaxies (e.g., iLisker et al.l 120061: 
Michielsen et al.l [2008 iChilingarianI 120091 : iSpolaor et al l 


2010al iRvs et al.l 201 .'^). as w e ll as i n many more massive 
ETGs bv iKuntschner et al.l (j20inl l. This finding adds 
further support to the notion of a continuity among low- 
mass and high-mass ETGs. 

5.4. Formation scenarios 

Over the past decade, many studies have argued for 
a scenario in which the low-mass (“dwarf’) ETGs in 
clusters originated from late-type, star-forming galax- 
ies that were tran s formed by the clu s ter environment 


(IMoore et alJ 1998 ; iBoselli et al.l 120081 IRvs et abi 12014 
iToloba et all 12001 12014^ Two main mechanisms 
are usually invoked to transform the infalling objects. 
First, ram-pressure stripping removes gas in about 100- 
200 Myr and rapidly quenches the star formation on a 
time-scale of ~ 1 Gyr (IBoselli et al.ll200^ . Second, grav¬ 
itational tidal harassment dynamically heats the low- 
mass ETG progenitors, lowering their angular momen¬ 
tum and concentrating the stellar comp onent (leading to 
more compact systems; IRvs et 301201411 . A third proces s 
is certainly also at play: starvation (ILarson et al.lll980ll . 
which is a form of gas stripping b ut acts on a much lo nger 
time scale (>5 Gyr) as shown bv IBoselli et al.ll2009l (see 
their Fig. 4), and affects more efficiently the gaseous halo 
rather than the inner regions. This process tends to sig¬ 
nificantly decrease the surface brightness of the affected 
object but keeps its characteristic size rather unchanged, 
in contra st with the effect o f ram-pressure stripping (see 
Fig. 3 of IBoselli et al.ll20(i^ . Starvation is therefore not 
able to solely explain the properties of low-mass ETGs, 
especially for the more compact ones. However it cer¬ 
tainly helps stopping their star formation, especially in 
clusters where a high rate of interactions can easily re¬ 


move the gas ([Rasmussen et al.ll2012ll . 

The exact nature of the prog enitors is still a matter 
of debate. ILisker et al.l (120131 1. based on a compari¬ 
son between Virgo observations and Millenium B sim¬ 
ulations, pointed out that late-type galaxies in the local 
universe may not be representative of the progenitors of 
today’s low-mass ETGs. They computed the time spent 
by a low-mass ETG inside a high-mass {M > 10^^ M 0 ) 
halo, as well as the mass-loss since first entering the halo 
(see their Fig. 4). According to their analysis, and us¬ 
ing the gala xy d ensity parameter Eio derived from the 
NGVS (see ? I5.1|) . our sample galaxies would have spent, 
on average, ^ 7 Gyr inside a halo more massive than 
10^^ Mq . Since our objects have a median age of 9 Gyr, 
this suggests that they have spent most of their lives 
inside a massive halo, experiencing a variety of environ¬ 
mental effects (i.e. from close pair to group interactions) 
that would certainly have helped shape their present day 
structure. Because the scatter in these correlations is 
quite large, it is hardly surprising for some galaxies to 
undergo infall at later times and/or experience quite dif¬ 
ferent interaction histories. This stochasticity could ex¬ 
plain the diversity in kinematics and stellar populations 
observed in low-mass ETGs, as well as the highly com¬ 
pact nature of some (rare) objects in dense regions. 

The galaxies in our sample have Xr values that are 
lower than more massive galaxies but consistent with 
“typical” values found for low-mass ETGs. Gravitational 
harassment (including tidal stripping) could explain such 
a decrease in angular momentum and one might naively 
expect that the longer a given galaxy stays in the clus¬ 
ter, the more angular momentum it looses and the more 
compact it becomes. However, our galaxy sample does 
not show a systematically lower specific angular momen¬ 
tum than less compact, low-mass ETGs. This may sug¬ 
gest either that their compactness is not related solely to 
gravitational harassment, or that this process (especially 
tidal stripping) can sometimes raise the stellar angular 
momentum of the system. 

Several of our program galaxies harbor younger and 
more metal-rich stellar populations in their center com¬ 
pared to their outskirts. This is consistent with what 
one would expect from ram-pressure stripping, where 
the pressure of the intra-cluster medium can remove gas 
from a galaxy but its efficiency depends on the relative 
strength of the gravitational potential of the object. The 
gas is thus more easily removed from the outskirts of a 
galaxy than from its center where it ca n stay bound to 
its ho st and continue forming stars fsee iFumagalli et al.l 
12014 in the case of a late-type galaxy). Gravitational 
tidal forces could also be responsible for younger stel¬ 
lar populations in the centers of low-mass ETGs as it 
promotes the formation of bars (iMaver et al.ll200^ that 
are known to induce gas infall. Others scenarios could 
potentially explain these observations, though: simple 
secular evolution, external accretion of gas, i.e. from the 
intergalactic medium or from a gas-rich companion (in¬ 
teraction, merger). The latter is more likely to occur in 
small groups, before a galaxy enters the dense regions of 
the present-day cluster where the relative velocities are 
high. 

The existence of 2 — cr galaxies in the low-mass ETGs 
population may also provide some constraints pertaining 
to their formation and evolution. The counter-rotating 
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disk in VCC 1475, as revealed by our two-dimensional 
spectroscopy, appears as a rather thin, central structure. 
N-body simulations predict that ^ 1% of the dwarfs 
in Vir go could exhibit counter rot ation due to harass¬ 
ment (|Gonzalez-Garcfa et al.l 120051) . However, the for¬ 
mation of such a counter-rotating structure via this pro¬ 
cess would make the outer part of the galaxy (i? > Re) 
counter rotating, and this is not what is observed for 
VCC 1475. A pure stellar merger event is also unlikely 
in this case as the inner disk is more metal-rich than the 
surrounding stellar component. A gas-accretion event 
is a more promising hypothesis, probably the result of 
an int eraction with a gas-rich galaxy (see iToloba et al.l 
l2014ai for more details about KDC formation scenarios 
in cluster environments). The exact origin of this 2 — cr 
galaxy is unclear but its appearance in a low-mass ETC 
sample certainly adds a level of complexity to the strip¬ 
ping and harassment scenario, emphasizing the impor¬ 
tance of the pre-cluster environmental history of cluster 
galaxies. 

In summary, we have shown that the properties of 
the stellar populations in low-mass ETCs show a level 
of continuity with more massive ETCs drawn from the 
ATLAS^^ survey (which itself is representative of the 
local volume within 40 Mpc). The age and metallicity 
follow the lines of iso-velocity dispersion (with a large 
scatter) suggesting that the formation and evolution pro¬ 
cesses that makes low-mass ETCs (whether compact or 
not) are not very different from the ones acting on high- 
mass ETCs, or at least that their impact vary smoothly 
with the mass, size and environment of the galaxy. 

6 . SUMMARY 

We have presented and analyzed spectroscopic obser¬ 
vations for eight low-mass compact ETCs in the Virgo 
cluster using the GMOS-IFU instrument mounted on 
the Gemini North 8 -m teles c ope. These objects were 
selected from iBinggeli et al.l (I1985D and our own Next 
Generation Virgo Survey fNGVS: Iherrarese et al.l[ 2 (H^ 
to probe a regime of ETG size and surface brightness 
that has largely been unexplored with IFU instruments. 
Our goal has been to characterize their basic kinematics 
and stellar populations in order to examine possible con¬ 
nections between this rare class of galaxy with more mas¬ 
sive (M > 10^° M©) and less compact low-mass ETGs. 
We performed kinematics and stellar populations anal¬ 
ysis using full spectral fitting. Our main results can be 
summarized as follows: 

• We calculated the dynamical mass inside one effec¬ 
tive radius for each of our objects and confirmed 
that these galaxies are more compact than most 
ETGs of comp arable mass (|Toloba et al.l l2014d : 
iRvs et aO 120141 1. The sample galaxies span a dy¬ 
namical mass range of 4.5 x 10® to 1.1 x 10^® Mq 
and have effective radii of 165 to 740 pc. Using the 
parameters L^ax/Lgai and E 15 as proxies of the 
local galaxy density, we have found that the most 
compact low-mass ETGs are preferentially found 
in high density regions and often close to a very 
massive galaxy. 

• The two-dimensional maps of stellar radial velocity, 
y, and velocity dispersion, a, of our objects show 


a large variety of kinematic structures: from non¬ 
rotating systems to those with well organized and 
high rotation, as well as velocity dispersion profiles 
that vary from flat to centrally peaked. Our data 
also provide unambiguous evidence for the 2 — cr 
nature of the low-mass galaxy VGG 1475. In most 
cases, low-mass compact ETGs with organized and 
high radial stellar velocity exhibit disky isophotes 
as revealed by the analysis of the NGVS deep im¬ 
ages. 

• We assessed the rotational support in these galax¬ 
ies by calculating their specific apparent angu¬ 
lar momentum, Aij^/ 2 ) finding two slow-rotators 
and six fast-rotators. Our sample galaxies have 
Afle /2 values similar to those found previously in 
more extended low -mass ETGs (jRvs et al.l I 201 .H 
iToloba et al.ll20143l and none of them reach values 
higher than 0.4 a s is sometimes observ ed in more 
massive galaxies (lEmsellem et al.ll2011[l . They are 
generally round objects and we do not find very 
flattened galaxies (e > 0.4). Three of the five most 
compact low-mass ETGs have disky isophotes, and 
all of them are fast rotators with the highest spe¬ 
cific angular momentum of our sample. 

• We found that our program galaxies contain old 
stars, between 6 and 11.5 Gyr, with metallicity be¬ 
tween —0.7 and 0.2. Half of our sample galaxies 
contain younger and more metal-rich stellar pop¬ 
ulations in their centers. Interestingly, low-mass 
compact ETGs with underlying disky isophotes are 
found to have the oldest stellar populations among 
our sample. 

• We derived SSP equivalent ages and metallicities 
within Re/2 for our program galaxies, finding a 
striking conti nuity between the low-m ass and high- 
mass ETGs (iGaopellari et al.l l2013al l. The more 
compact galaxies, at a fixed mass, are found to 
have more metal-rich and older (albeit with a larger 
scatter) stell ar populations. We also extended the 
results from iMcDermid et al.l (I2015D to low-mass 
ETGs, showing that the velocity dispersion is a 
relatively good predictor of the stellar population 
properties. 

Our results are consistent with a scenario in which 
low-mass, compact ETGs have spent a large fraction 
of their lives in the cluster environment and have been 
strongly shaped by environmental mechanisms such as 
ram-pressure stripping, tidal stripping, starvation and 
gravitational harassment. This picture is quantitatively 
supported b y a com parison of our data to the study of 
iLisker et al.l (I2013I1 . who jointly analyzed Millenium H 
simulations and Virgo galaxy observations. Our study 
underscores the importance of environmental history for 
low-mass ETGs, such as close-pairs, groups, etc. before 
the pre-present day cluster epoch. Perhaps most signifi¬ 
cantly, our study adds to a growing body of evidence that 
a continuity exists in the structural, kinematic and stel¬ 
lar population properties of low- and high-mass ETGs, 
which in turn suggests that ETG formation and evolu¬ 
tion processes vary smoothly with mass, size and envi- 
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ronment across the full range in galaxy mass. A more de¬ 
tailed account of the relative weight of each evolutionary 
processes should wait for observations of larger galaxy 
samples and realistic numerical simulations. 

7. ACKNOWLEDGEMENTS 

This work is based on observations obtained at the 
Gemini Observatory processed using the Gemini IRAF 
package, which is operated by the Association of Univer¬ 
sities for Research in Astronomy, Inc., under a cooper¬ 
ative agreement with the NSF on behalf of the Gemini 
partnership: the National Science Foundation (United 
States), the National Research Gouncil (Ganada), CON- 
ICYT (Ghile), the Australian Research Gouncil (Aus¬ 
tralia), Ministrio da Cincia, Tecnologia e Inovao (Brazil) 
and Ministerio de Giencia, Tecnologa e Innovacin Pro- 
ductiva (Argentina). 

The authors thank Elisa Toloba for sharing her work 
ahead of publication and Mina Koleva for the use of her 
data for VCC 1475. 

AG acknowledges the hospitality and support of the 


Allington-Smith, J., Murray, G., Content, R., et al. 2002, Publ. 
Astron. Soc. Pacific, 114, 892 

Arrigoni Battaia, F., Gavazzi, G., Fumagalli, M., et al. 2012, A&A, 
543, A112 

Bacon, R., Copin, Y., Monnet, G., et al. 2001, MNRAS, 326, 23 
Bender, R., Burstein, D., & Faber, S. M. 1992, ApJ, 399, 462 
Berlin, G., Ciotti, L., &c Del Principe, M. 2002, A&A, 386, 149 
Binggeli, B., Sandage, A., & Tammann, G. A. 1985, AJ, 90, 1681 
Blakeslee, J. P., Jordan, A., Mei, S., et al. 2009, ApJ, 694, 556 
Boselli, A., Boissier, S., Cortese, L., & Gavazzi, G. 2008, A&A, 
489, 1015 

Boselli, a., Boissier, S., Cortese, L., & Gavazzi, G. 2009, Astron. 
Nachrichten, 330, 904 

Boselli, A., Voyer, E., Boissier, S., et al. 2014, arXiv:1407.4986 
Cappellari, M., & Copin, Y. 2003, MNRAS, 342, 345 
Cappellari, M., & Emsellem, E. 2004, Publ. Astron. Soc. Pacific, 
116, 138 

Cappellari, M., Bacon, R., Bureau, M., et al. 2006, MNRAS, 366, 
1126 

Cappellari, M., Emsellem, E., Krajnovic, D., et al. 2011, MNRAS, 

413, 813 

Cappellari, M., Scott, N., Alatalo, K., et al. 2013a, MNRAS, 432, 
1709 

Cappellari, M., McDermid, R. M., Alatalo, K., et al. 2013b, 
MNRAS, 432, 1862 

Cenarro, A. J., Cardiel, N., Gorgas, J., et al. 2001, MNRAS, 326, 
959 

Chilingarian, I., Cayatte, V., Revaz, Y., et al. 2009, Science, 326, 
1379 

Chilingarian, I. V. 2009, MNRAS, 394, 1229 

Cote, P., Blakeslee, J. P., Ferrarese, L., et al. 2004, ApJS, 153, 223 
Cote, P., Piatek, S., Ferrarese, L., et al. 2006, ApJS, 165, 57 
Cote, P., Ferrarese, L., Jordan, A., et al. 2007, ApJ, 671, 1456 
Dabringhausen, J., Hilker, M., & Kroupa, P. 2008, MNRAS, 386, 
864 

de Zeeuw, P., Bureau, M., Emsellem, E., et al. 2002, MNRAS, 329, 
513 

Durrell, P. R., Cote, P., Peng, E. W., et al. 2014, ApJ, 794, 103 
Emsellem, E., Cappellari, M., Krajnovi, D., et al. 2007, MNRAS, 
379, 401 

Emsellem, E., Cappellari, M., Krajnovic, D., et al. 2011, MNRAS, 

414, 888 

Faber, S. M. 1973, ApJ, 179, 423 

Falcon-Barroso, J., Sanchez-Blazquez, P., Vazdekis, A., et al. 2011, 
A&A, 532, A95 

Ferrarese, L., Cote, P., Dalla Bonta, E., et al. 2006, ApJ, 644, L21 
Ferrarese, L., Cote, P., Cuillandre, J.-C., et al. 2012, ApJS, 200, 4 
Fumagalli, M., Fossati, M., Hau, G. K. T., et al. 2014, 9, 9 
Gavazzi, G., Donati, a., Cucciati, O., et al. 2005, A&A, 430, 411 


Macquarie University Astronomy and Astrophysics Gen- 
ter during the preparation of this work. 

EWP acknowledges support from the National Natural 
Science Foundation of China under Grant No. 11173003, 
and from the Strategic Priority Research Program, “The 
Emergence of Cosmological Structures”, of the Chinese 
Academy of Sciences, Grant No. XDB09000105. 

This work is supported in part by the French Agence 
Nationale de la Recherche (ANR) Grant Programme 
Blanc VIRAGE (ANRlO-BLANC-0506-01), and by the 
Canadian Advanced Network for Astronomical Research 
(CANFAR) which has been made possible by funding 
from CANARIE under the Network-Enabled Platforms 
program. This research used the facilities of the Cana¬ 
dian Astronomy Data Center operated by the National 
Research Council of Canada with the support of the 
Canadian Space Agency. The NGVS team owes a debt 
of gratitude to the director and the staff of the Canada- 
France-Hawaii Telescope, whose dedication, ingenuity, 
and expertise have helped make the survey a reality. 

Facilities: Gemini, CFHT 


Gnedin, O. Y. 2003, ApJ, 582, 141 

Gonzalez-Garci'a, a. G., Aguerri, J. a. L., Balcells, M. 2005, A&A, 
444, 803 

Graham, A. W. 2002, ApJ, 568, L13 
Graham, A. W., & Guzman, R. 2003, AJ, 125, 2936 
Graves, G. J., & Faber, S. M. 2010, ApJ, 717, 803 
Hook, 1. M., J0 rgensen, L, AllingtonSmith, J. R., et al. 2004, Publ. 
Astron. Soc. Pacific, 116, 425 

Howley, K. M., Guhathakurta, P., van der Marel, R., et al. 2013, 
ApJ, 765, 65 

Huxor, a. P., Phillipps, S., Price, J., & Harniman, R. 2011, 
MNRAS, 414, 3557 

Jedrzejewski, R. I. 1987, in lAU Symposium, Vol. 127, Structure 
and Dynamics of Elliptical Galaxies, ed. P. T. de Zeeuw, 37—44 
Jerjen, H., & Binggeli, B. 1997, in Astronomical Society of the 
Pacific Gonference Series, Vol. 116, The Nature of Elliptical 
Galaxies; 2nd Stromlo Symposium, ed. M. Arnaboldi, G. S. Da 
Costa, & P. Saha, 239 

Jordan, A., McLaughlin, D. E., Cote, P., et al. 2007, ApJS, 171, 
101 

Keenan, D. W., & Innanen, K. A. 1975, AJ, 80, 290 
King, I. 1962, AJ, 67, 274 

Koleva, M., Bedregal, A., Prugniel, P., De Rijcke, S., & Zeilinger, 
W. 2011a, EAS Publ. Ser., 48, 85 
Koleva, M., Prugniel, P., De Rijcke, S., & Zeilinger, W. W. 2011b, 
MNRAS, 417, 1643 

Koleva, M., Prugniel, P., Ocvirk, P., Le Borgne, D., & Soubiran, 
C. 2008, MNRAS, 385, 1998 
Kormendy, J. 1985, ApJ, 295, 73 
Kormendy, J., & Bender, R. 1997, ApJ, 3 

Kormendy, J., Fisher, D. B., Cornell, M. E., & Bender, R. 2009, 
ApJS, 182, 216 

Krajnovic, D., Emsellem, E., Cappellari, M., et al. 2011, MNRAS, 
414, 2923 

Kuntschner, H., Emsellem, E., Bacon, R., et al. 2010, MNRAS, 
408, 97 

Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 
692 

Lauer, T. R., Tremaine, S., Ajhar, E. A., et al. 1996, ApJ, 471, L79 
Lisker, T., Glatt, K., Westera, P., & Grebel, E. K. 2006, AJ, 132, 
2432 

Lisker, T., Grebel, E. K., Binggeli, B., & Glatt, K. 2007, ApJ, 660, 
1186 

Lisker, T., Weinmann, S. M., Janz, J., & Meyer, H. T. 2013, 
MNRAS, 432, 1162 

Mayer, L., Mastropietro, C., Wadsley, J., Stadel, J., & Moore, B. 
2006, MNRAS, 369, 1021 

McDermid, R. M., Alatalo, K., Blitz, L., et al. 2015, ArXiv e-prints, 
arXiv:1501.03723 


NGVS. XII. Stellar Populations and Kinematics of Compact, Low-Mass Early-Type Galaxies 


19 


Mei, S., Blakeslee, J. P., Cote, P., et al. 2007, ApJ, 655, 144 
Metz, M., & Kroupa, P. 2007, MNRAS, 376, 387 
Michielsen, D., Koleva, M., Prugniel, P., et al. 2007, ApJ, 670, LlOl 
Michielsen, D., Boselli, a., Conselice, C. J., et al. 2008, MNRAS, 
385, 1374 

Monachesi, A., Trager, S. C., Lauer, T. R., et al. 2012, ApJ, 745, 
97 

Moore, B., Lake, G., Katz, N. 1998, ApJ, 495, 139 

Munoz, R. P., Puzia, T. H., Langon, A., et al. 2014, ApJS, 210, 4 

Okazaki, T., &; Taniguchi, Y. 2000, 1, 9 

Paudel, S., Lisker, T., & Kuntschner, H. 2011, MNRAS, 413, 1764 
Paudel, S., Due, P.-A., Cote, P., et al. 2013, ApJ, 767, 133 
Peng, E. W., Jordan, A., Cote, P., et al. 2008, ApJ, 681, 197 
Prugniel, P., Chilingarian, I., Sil’chenko, O., Sz Afanasiev, V. 2005, 
Proc. Int. Astron. Union, 1, 73 

Rasmussen, J., Mulchaey, J. S., Bai, L., et al. 2012, ApJ, 757, 122 
Rys, A., Falcon-Barroso, J., Sz van de Yen, G. 2013, MNRAS, 428, 
2980 

Rys, A., van de Ven, G., &; Falcon-Barroso, J. 2014, MNRAS, 439, 
284 

Salpeter, E. E. 1955, ApJ, 121, 161 

Sanchez, S. F., Kennicutt, R. C., Gil de Paz, A., et al. 2012, A&A, 
538, A8 

Sanchez-Blazquez, P., Gorgas, J., Cardiel, N., Sz Gonzalez, J. J. 
2006, A&A, 457, 809 

Sanchez-Blazquez, P., Peletier, R. F., Jimenez-Vicente, J., et al. 
2006, MNRAS, 371, 703 


Seth, A. C., van den Bosch, R., Mieske, S., et al. 2014, Natur, 513, 
398 

Smith, R. J., Lucey, J. R., Price, J., Hudson, M. J., Sz Phillipps, S. 
2012, MNRAS, 419, 3167 

Spolaor, M., Hau, G. K. T., Forbes, D. a., Sz Couch, W. J. 2010a, 
MNRAS, 408, 254 

Spolaor, M., Kobayashi, C., Forbes, D. a.. Couch, W. J., Sz Hau, 
G. K. T. 2010b, MNRAS, 408, 272 
Thomas, D., Maraston, C., Schawinski, K., Sarzi, M., Sz Silk, J. 
2010, MNRAS, 1789, 1775 

Toloba, E., Boselli, A., Cenarro, A. J., et al. 2011, A&A, 526, A114 
Toloba, E., Boselli, a., Gorgas, J., et al. 2009, ApJ, 707, L17 
Toloba, E., Guhathakurta, P., van de Ven, G., et al. 2014a, ApJ, 
783, 120 

Toloba, E., Guhathakurta, P., Peletier, R., et al. 2014b, 

arXiv: 1410.1550 

Toloba, E., Guhathakurta, P., Boselli, A., et al. 2014c, 

arXiv: 1410.1552 

Turner, M. L., Gote, P., Ferrarese, L., et al. 2012, ApJS, 203, 5 
Vazdekis, a., Ricciardelli, E., Cenarro, a. J., et al. 2012, MNRAS, 
424, 157 

Villegas, D., Jordan, A., Peng, E. W., et al. 2010, ApJ, 717, 603 
Vollmer, B., Soida, M., Chung, A., et al. 2009, 143, 8 
Wake, D. a., van Dokkum, P. G., & Franx, M. 2012, ApJ, 751, L44 
Zhu, L., Long, R. J., Mao, S., et al. 2014, ApJ, 792, 59 


